TLRs: Innate Immune Sentries against SARS-CoV-2 Infection.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been responsible for a devastating pandemic since March 2020. Toll-like receptors (TLRs), crucial components in the initiation of innate immune responses to different pathogens, trigger the downstream production of pro-inflammatory cytokines, interferons, and other mediators. It has been demonstrated that they contribute to the dysregulated immune response observed in patients with severe COVID-19. TLR2, TLR3, TLR4 and TLR7 have been associated with COVID-19 severity. Here, we review the role of TLRs in the etiology and pathogenesis of COVID-19, including TLR7 and TLR3 rare variants, the L412F polymorphism in TLR3 that negatively regulates anti-SARS-CoV-2 immune responses, the TLR3-related cellular senescence, the interaction of TLR2 and TLR4 with SARS-CoV-2 proteins and implication of TLR2 in NET formation by SARS-CoV-2. The activation of TLRs contributes to viral clearance and disease resolution. However, TLRs may represent a double-edged sword which may elicit dysregulated immune signaling, leading to the production of proinflammatory mediators, resulting in severe disease. TLR-dependent excessive inflammation and TLR-dependent antiviral response may tip the balance towards the former or the latter, altering the equilibrium that drives the severity of disease.

SELP Asp603Asn and severe thrombosis in COVID-19 males.

Thromboembolism is a frequent cause of severity and mortality in COVID-19. However, the etiology of this phenomenon is not well understood. A cohort of 1186 subjects, from the GEN-COVID consortium, infected by SARS-CoV-2 with different severity was stratified by sex and adjusted by age. Then, common coding variants from whole exome sequencing were mined by LASSO logistic regression. The homozygosity of the cell adhesion molecule P-selectin gene (SELP) rs6127 (c.1807G > A; p.Asp603Asn) which has been already associated with thrombotic risk is found to be associated with severity in the male subcohort of 513 subjects (odds ratio = 2.27, 95% Confidence Interval 1.54-3.36). As the SELP gene is downregulated by testosterone, the odd ratio is increased in males older than 50 (OR 2.42, 95% CI 1.53-3.82). Asn/Asn homozygotes have increased D-dimers values especially when associated with poly Q ≥ 23 in the androgen receptor (OR 3.26, 95% CI 1.41-7.52). These results provide a rationale for the repurposing of antibodies against P-selectin as adjuvant therapy in rs6127 male homozygotes especially if older than 50 or with an impaired androgen receptor.

Gain- and Loss-of-Function CFTR Alleles Are Associated with COVID-19 Clinical Outcomes.

Carriers of single pathogenic variants of the CFTR (cystic fibrosis transmembrane conductance regulator) gene have a higher risk of severe COVID-19 and 14-day death. The machine learning post-Mendelian model pinpointed CFTR as a bidirectional modulator of COVID-19 outcomes. Here, we demonstrate that the rare complex allele [G576V;R668C] is associated with a milder disease via a gain-of-function mechanism. Conversely, CFTR ultra-rare alleles with reduced function are associated with disease severity either alone (dominant disorder) or with another hypomorphic allele in the second chromosome (recessive disorder) with a global residual CFTR activity between 50 to 91%. Furthermore, we characterized novel CFTR complex alleles, including [A238V;F508del], [R74W;D1270N;V201M], [I1027T;F508del], [I506V;D1168G], and simple alleles, including R347C, F1052V, Y625N, I328V, K68E, A309D, A252T, G542*, V562I, R1066H, I506V, I807M, which lead to a reduced CFTR function and thus, to more severe COVID-19. In conclusion, CFTR genetic analysis is an important tool in identifying patients at risk of severe COVID-19.

An explainable model of host genetic interactions linked to COVID-19 severity.

We employed a multifaceted computational strategy to identify the genetic factors contributing to increased risk of severe COVID-19 infection from a Whole Exome Sequencing (WES) dataset of a cohort of 2000 Italian patients. We coupled a stratified k-fold screening, to rank variants more associated with severity, with the training of multiple supervised classifiers, to predict severity based on screened features. Feature importance analysis from tree-based models allowed us to identify 16 variants with the highest support which, together with age and gender covariates, were found to be most predictive of COVID-19 severity. When tested on a follow-up cohort, our ensemble of models predicted severity with high accuracy (ACC = 81.88%; AUCROC = 96%; MCC = 61.55%). Our model recapitulated a vast literature of emerging molecular mechanisms and genetic factors linked to COVID-19 response and extends previous landmark Genome-Wide Association Studies (GWAS). It revealed a network of interplaying genetic signatures converging on established immune system and inflammatory processes linked to viral infection response. It also identified additional processes cross-talking with immune pathways, such as GPCR signaling, which might offer additional opportunities for therapeutic intervention and patient stratification. Publicly available PheWAS datasets revealed that several variants were significantly associated with phenotypic traits such as "Respiratory or thoracic disease", supporting their link with COVID-19 severity outcome.

Multiomic analysis reveals cell-type-specific molecular determinants of COVID-19 severity.

The determinants of severe COVID-19 in healthy adults are poorly understood, which limits the opportunity for early intervention. We present a multiomic analysis using machine learning to characterize the genomic basis of COVID-19 severity. We use single-cell multiome profiling of human lungs to link genetic signals to cell-type-specific functions. We discover >1,000 risk genes across 19 cell types, which account for 77% of the SNP-based heritability for severe disease. Genetic risk is particularly focused within natural killer (NK) cells and T cells, placing the dysfunction of these cells upstream of severe disease. Mendelian randomization and single-cell profiling of human NK cells support the role of NK cells and further localize genetic risk to CD56bright NK cells, which are key cytokine producers during the innate immune response. Rare variant analysis confirms the enrichment of severe-disease-associated genetic variation within NK-cell risk genes. Our study provides insights into the pathogenesis of severe COVID-19 with potential therapeutic targets.

Pathogen-sugar interactions revealed by universal saturation transfer analysis.

Many pathogens exploit host cell-surface glycans. However, precise analyses of glycan ligands binding with heavily modified pathogen proteins can be confounded by overlapping sugar signals and/or compounded with known experimental constraints. Universal saturation transfer analysis (uSTA) builds on existing nuclear magnetic resonance spectroscopy to provide an automated workflow for quantitating protein-ligand interactions. uSTA reveals that early-pandemic, B-origin-lineage severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike trimer binds sialoside sugars in an "end-on" manner. uSTA-guided modeling and a high-resolution cryo-electron microscopy structure implicate the spike N-terminal domain (NTD) and confirm end-on binding. This finding rationalizes the effect of NTD mutations that abolish sugar binding in SARS-CoV-2 variants of concern. Together with genetic variance analyses in early pandemic patient cohorts, this binding implicates a sialylated polylactosamine motif found on tetraantennary N-linked glycoproteins deep in the human lung as potentially relevant to virulence and/or zoonosis.

Geographical distribution of cystic fibrosis carriers as population genetic determinant of COVID-19 spread and fatality in 37 countries.

COVID-19 has shown a relevant heterogeneity in spread and fatality among countries together with a significant variability in its clinical presentation, indicating that host genetic factors may influence COVID-19 pathogenicity. Indeed, subjects carrying single pathogenic variants of the Cystic Fibrosis (CF) Transmembrane Conductance Regulator (CFTR) gene - i.e. CF carriers - are more susceptible to respiratory tract infections and are more likely to undergo severe COVID-19 with higher risk of 14-day mortality. Given that CF carrier prevalence varies among ethnicities and nations, an ecological study in 37 countries was conducted, in order to determine to what extent the diverse CF carrier geographical distribution may have affected COVID-19 spread and fatality during the first pandemic wave. The CF prevalence in countries, as indicator of the geographical distribution of CF carriers, significantly correlated in a direct manner with both COVID-19 prevalence and its Case Fatality Rate (CFR). In a regression study weighted for the number of tests performed, COVID-19 prevalence positively correlated with CF prevalence, while CFR correlated with population percentage older than 65-year, cancer and CF prevalence. Multivariate regression model also confirmed COVID-19 CFR to be associated with CF prevalence, after adjusting for elderly, cancer prevalence, and weighting for the number of tests performed. This study suggests a putative contribution of population genetics of CFTR in understanding the spatial distribution of COVID-19 spread and fatality.

Carriers of ADAMTS13 Rare Variants Are at High Risk of Life-Threatening COVID-19.

Thrombosis of small and large vessels is reported as a key player in COVID-19 severity. However, host genetic determinants of this susceptibility are still unclear. Congenital Thrombotic Thrombocytopenic Purpura is a severe autosomal recessive disorder characterized by uncleaved ultra-large vWF and thrombotic microangiopathy, frequently triggered by infections. Carriers are reported to be asymptomatic. Exome analysis of about 3000 SARS-CoV-2 infected subjects of different severities, belonging to the GEN-COVID cohort, revealed the specific role of vWF cleaving enzyme ADAMTS13 (A disintegrin-like and metalloprotease with thrombospondin type 1 motif, 13). We report here that ultra-rare variants in a heterozygous state lead to a rare form of COVID-19 characterized by hyper-inflammation signs, which segregates in families as an autosomal dominant disorder conditioned by SARS-CoV-2 infection, sex, and age. This has clinical relevance due to the availability of drugs such as Caplacizumab, which inhibits vWF-platelet interaction, and Crizanlizumab, which, by inhibiting P-selectin binding to its ligands, prevents leukocyte recruitment and platelet aggregation at the site of vascular damage.

Host genetic basis of COVID-19: from methodologies to genes.

The COVID-19 pandemic caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is having a massive impact on public health, societies, and economies worldwide. Despite the ongoing vaccination program, treating COVID-19 remains a high priority; thus, a better understanding of the disease is urgently needed. Initially, susceptibility was associated with age, sex, and other prior existing comorbidities. However, as these conditions alone could not explain the highly variable clinical manifestations of SARS-CoV-2 infection, the attention was shifted toward the identification of the genetic basis of COVID-19. Thanks to international collaborations like The COVID-19 Host Genetics Initiative, it became possible the elucidation of numerous genetic markers that are not only likely to help in explaining the varied clinical outcomes of COVID-19 patients but can also guide the development of novel diagnostics and therapeutics. Within this framework, this review delineates GWAS and Burden test as traditional methodologies employed so far for the discovery of the human genetic basis of COVID-19, with particular attention to recently emerged predictive models such as the post-Mendelian model. A summary table with the main genome-wide significant genomic loci is provided. Besides, various common and rare variants identified in genes like TLR7, CFTR, ACE2, TMPRSS2, TLR3, and SELP are further described in detail to illustrate their association with disease severity.

The polymorphism L412F in TLR3 inhibits autophagy and is a marker of severe COVID-19 in males.

The polymorphism L412F in TLR3 has been associated with several infectious diseases. However, the mechanism underlying this association is still unexplored. Here, we show that the L412F polymorphism in TLR3 is a marker of severity in COVID-19. This association increases in the sub-cohort of males. Impaired macroautophagy/autophagy and reduced TNF/TNFα production was demonstrated in HEK293 cells transfected with TLR3L412F-encoding plasmid and stimulated with specific agonist poly(I:C). A statistically significant reduced survival at 28 days was shown in L412F COVID-19 patients treated with the autophagy-inhibitor hydroxychloroquine (p = 0.038). An increased frequency of autoimmune disorders such as co-morbidity was found in L412F COVID-19 males with specific class II HLA haplotypes prone to autoantigen presentation. Our analyses indicate that L412F polymorphism makes males at risk of severe COVID-19 and provides a rationale for reinterpreting clinical trials considering autophagy pathways.Abbreviations: AP: autophagosome; AUC: area under the curve; BafA1: bafilomycin A1; COVID-19: coronavirus disease-2019; HCQ: hydroxychloroquine; RAP: rapamycin; ROC: receiver operating characteristic; SARS-CoV-2: severe acute respiratory syndrome coronavirus 2; TLR: toll like receptor; TNF/TNF-α: tumor necrosis factor.

Rare variants in Toll-like receptor 7 results in functional impairment and downregulation of cytokine-mediated signaling in COVID-19 patients.

Toll-like receptors (TLR) are crucial components in the initiation of innate immune responses to a variety of pathogens, triggering the production of pro-inflammatory cytokines and type I and II interferons, which are responsible for innate antiviral responses. Among the different TLRs, TLR7 recognizes several single-stranded RNA viruses including SARS-CoV-2. We and others identified rare loss-of-function variants in X-chromosomal TLR7 in young men with severe COVID-19 and with no prior history of major chronic diseases, that were associated with impaired TLR7 signaling as well as type I and II IFN responses. Here, we performed RNA sequencing to investigate transcriptome variations following imiquimod stimulation of peripheral blood mononuclear cells isolated from patients carrying previously identified hypomorphic, hypofunctional, and loss-of-function TLR7 variants. Our investigation revealed a profound impairment of the TLR7 pathway in patients carrying loss-of-function variants. Of note, a failure in IFNγ upregulation following stimulation was also observed in cells harboring the hypofunctional and hypomorphic variants. We also identified new TLR7 variants in severely affected male patients for which a functional characterization of the TLR7 pathway was performed demonstrating a decrease in mRNA levels in the IFNα, IFNγ, RSAD2, ACOD1, IFIT2, and CXCL10 genes.

Common, low-frequency, rare, and ultra-rare coding variants contribute to COVID-19 severity.

The combined impact of common and rare exonic variants in COVID-19 host genetics is currently insufficiently understood. Here, common and rare variants from whole-exome sequencing data of about 4000 SARS-CoV-2-positive individuals were used to define an interpretable machine-learning model for predicting COVID-19 severity. First, variants were converted into separate sets of Boolean features, depending on the absence or the presence of variants in each gene. An ensemble of LASSO logistic regression models was used to identify the most informative Boolean features with respect to the genetic bases of severity. The Boolean features selected by these logistic models were combined into an Integrated PolyGenic Score that offers a synthetic and interpretable index for describing the contribution of host genetics in COVID-19 severity, as demonstrated through testing in several independent cohorts. Selected features belong to ultra-rare, rare, low-frequency, and common variants, including those in linkage disequilibrium with known GWAS loci. Noteworthily, around one quarter of the selected genes are sex-specific. Pathway analysis of the selected genes associated with COVID-19 severity reflected the multi-organ nature of the disease. The proposed model might provide useful information for developing diagnostics and therapeutics, while also being able to guide bedside disease management.

Age-dependent impact of the major common genetic risk factor for COVID-19 on severity and mortality.

BackgroundThere is considerable variability in COVID-19 outcomes among younger adults, and some of this variation may be due to genetic predisposition.MethodsWe combined individual level data from 13,888 COVID-19 patients (n = 7185 hospitalized) from 17 cohorts in 9 countries to assess the association of the major common COVID-19 genetic risk factor (chromosome 3 locus tagged by rs10490770) with mortality, COVID-19-related complications, and laboratory values. We next performed metaanalyses using FinnGen and the Columbia University COVID-19 Biobank.ResultsWe found that rs10490770 risk allele carriers experienced an increased risk of all-cause mortality (HR, 1.4; 95% CI, 1.2-1.7). Risk allele carriers had increased odds of several COVID-19 complications: severe respiratory failure (OR, 2.1; 95% CI, 1.6-2.6), venous thromboembolism (OR, 1.7; 95% CI, 1.2-2.4), and hepatic injury (OR, 1.5; 95% CI, 1.2-2.0). Risk allele carriers age 60 years and younger had higher odds of death or severe respiratory failure (OR, 2.7; 95% CI, 1.8-3.9) compared with those of more than 60 years (OR, 1.5; 95% CI, 1.2-1.8; interaction, P = 0.038). Among individuals 60 years and younger who died or experienced severe respiratory failure, 32.3% were risk-variant carriers compared with 13.9% of those not experiencing these outcomes. This risk variant improved the prediction of death or severe respiratory failure similarly to, or better than, most established clinical risk factors.ConclusionsThe major common COVID-19 genetic risk factor is associated with increased risks of morbidity and mortality, which are more pronounced among individuals 60 years or younger. The effect was similar in magnitude and more common than most established clinical risk factors, suggesting potential implications for future clinical risk management.

C9orf72 Intermediate Repeats Confer Genetic Risk for Severe COVID-19 Pneumonia Independently of Age.

A cytokine storm, autoimmune features and dysfunctions of myeloid cells significantly contribute to severe coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Genetic background of the host seems to be partly responsible for severe phenotype and genes related to innate immune response seem critical host determinants. The C9orf72 gene has a role in vesicular trafficking, autophagy regulation and lysosome functions, is highly expressed in myeloid cells and is involved in immune functions, regulating the lysosomal degradation of mediators of innate immunity. A large non-coding hexanucleotide repeat expansion (HRE) in this gene is the main genetic cause of frontotemporal dementia (FTD) and amyotrophic lateral sclerosis (ALS), both characterized by neuroinflammation and high systemic levels of proinflammatory cytokines, while HREs of intermediate length, although rare, are more frequent in autoimmune disorders. C9orf72 full mutation results in haploinsufficiency and intermediate HREs seem to modulate gene expression as well and impair autophagy. Herein, we sought to explore whether intermediate HREs in C9orf72 may be a risk factor for severe COVID-19. Although we found intermediate HREs in only a small portion of 240 patients with severe COVID-19 pneumonia, the magnitude of risk for requiring non-invasive or mechanical ventilation conferred by harboring intermediate repeats >10 units in at least one C9orf72 allele was more than twice respect to having shorter expansions, when adjusted for age (odds ratio (OR) 2.36; 95% confidence interval (CI) 1.04-5.37, p = 0.040). The association between intermediate repeats >10 units and more severe clinical outcome (p = 0.025) was also validated in an independent cohort of 201 SARS-CoV-2 infected patients. These data suggest that C9orf72 HREs >10 units may influence the pathogenic process driving more severe COVID-19 phenotypes.

Severe COVID-19 in Hospitalized Carriers of Single CFTR Pathogenic Variants.

The clinical presentation of COVID-19 is extremely heterogeneous, ranging from asymptomatic to severely ill patients. Thus, host genetic factors may be involved in determining disease presentation and progression. Given that carriers of single cystic fibrosis (CF)-causing variants of the CFTR gene-CF-carriers-are more susceptible to respiratory tract infections, our aim was to determine their likelihood of undergoing severe COVID-19. We implemented a cohort study of 874 individuals diagnosed with COVID-19, during the first pandemic wave in Italy. Whole exome sequencing was performed and validated CF-causing variants were identified. Forty subjects (16 females and 24 males) were found to be CF-carriers. Among mechanically ventilated patients, CF-carriers were more represented (8.7%) and they were significantly (p < 0.05) younger (mean age 51 years) compared to noncarriers (mean age 61.42 years). Furthermore, in the whole cohort, the age of male CF-carriers was lower, compared to noncarriers (p < 0.05). CF-carriers had a relative risk of presenting an abnormal inflammatory response (CRP ≥ 20 mg/dL) of 1.69 (p < 0.05) and their hazard ratio of death at day 14 was 3.10 (p < 0.05) in a multivariate regression model, adjusted for age, sex and comorbidities. In conclusion, CF-carriers are more susceptible to the severe form of COVID-19, showing also higher risk of 14-day death.

Vitamin D and COVID-19 susceptibility and severity in the COVID-19 Host Genetics Initiative: A Mendelian randomization study.

BACKGROUND: Increased vitamin D levels, as reflected by 25-hydroxy vitamin D (25OHD) measurements, have been proposed to protect against COVID-19 based on in vitro, observational, and ecological studies. However, vitamin D levels are associated with many confounding variables, and thus associations described to date may not be causal. Vitamin D Mendelian randomization (MR) studies have provided results that are concordant with large-scale vitamin D randomized trials. Here, we used 2-sample MR to assess evidence supporting a causal effect of circulating 25OHD levels on COVID-19 susceptibility and severity. METHODS AND FINDINGS: Genetic variants strongly associated with 25OHD levels in a genome-wide association study (GWAS) of 443,734 participants of European ancestry (including 401,460 from the UK Biobank) were used as instrumental variables. GWASs of COVID-19 susceptibility, hospitalization, and severe disease from the COVID-19 Host Genetics Initiative were used as outcome GWASs. These included up to 14,134 individuals with COVID-19, and up to 1,284,876 without COVID-19, from up to 11 countries. SARS-CoV-2 positivity was determined by laboratory testing or medical chart review. Population controls without COVID-19 were also included in the control groups for all outcomes, including hospitalization and severe disease. Analyses were restricted to individuals of European descent when possible. Using inverse-weighted MR, genetically increased 25OHD levels by 1 standard deviation on the logarithmic scale had no significant association with COVID-19 susceptibility (odds ratio [OR] = 0.95; 95% CI 0.84, 1.08; p = 0.44), hospitalization (OR = 1.09; 95% CI: 0.89, 1.33; p = 0.41), and severe disease (OR = 0.97; 95% CI: 0.77, 1.22; p = 0.77). We used an additional 6 meta-analytic methods, as well as conducting sensitivity analyses after removal of variants at risk of horizontal pleiotropy, and obtained similar results. These results may be limited by weak instrument bias in some analyses. Further, our results do not apply to individuals with vitamin D deficiency. CONCLUSIONS: In this 2-sample MR study, we did not observe evidence to support an association between 25OHD levels and COVID-19 susceptibility, severity, or hospitalization. Hence, vitamin D supplementation as a means of protecting against worsened COVID-19 outcomes is not supported by genetic evidence. Other therapeutic or preventative avenues should be given higher priority for COVID-19 randomized controlled trials.

Association of Toll-like receptor 7 variants with life-threatening COVID-19 disease in males: findings from a nested case-control study.

Background: Recently, loss-of-function variants in TLR7 were identified in two families in which COVID-19 segregates like an X-linked recessive disorder environmentally conditioned by SARS-CoV-2. We investigated whether the two families represent the tip of the iceberg of a subset of COVID-19 male patients. Methods: This is a nested case-control study in which we compared male participants with extreme phenotype selected from the Italian GEN-COVID cohort of SARS-CoV-2-infected participants (<60 y, 79 severe cases versus 77 control cases). We applied the LASSO Logistic Regression analysis, considering only rare variants on young male subsets with extreme phenotype, picking up TLR7 as the most important susceptibility gene. Results: Overall, we found TLR7 deleterious variants in 2.1% of severely affected males and in none of the asymptomatic participants. The functional gene expression profile analysis demonstrated a reduction in TLR7-related gene expression in patients compared with controls demonstrating an impairment in type I and II IFN responses. Conclusions: Young males with TLR7 loss-of-function variants and severe COVID-19 represent a subset of male patients contributing to disease susceptibility in up to 2% of severe COVID-19. Funding: Funded by private donors for the Host Genetics Research Project, the Intesa San Paolo for 2020 charity fund, and the Host Genetics Initiative. Clinical trial number: NCT04549831.